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Abstract Triangle centers such as the center of gravity, the circumcenter, the orthocenter are well studied by 
geometers. Recently, under the guidance of Clark Kimberling, an electronic encyclopedia of triangle centers (ETC) 
has been developed, it contains more than 7000 centers and many properties of these points. In this paper, we 
describe how we created a certified version of ETC such that some of the properties described come along with a 
computer checked proof of its validity. 
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1 Introduction 


The characteristic points of triangles such as the center of gravity, the circumcenter, the orthocenter have been 
well studied for centuries by geometers. Some properties of these points are well known. For example the three 
points we cited are collinear, they belong to a line called the Euler line. More recently, under the guidance of Clark 
Kimberling [9-11], an electronic encyclopedia of triangle center (ETC) has been developed.' This encyclopedia 
contains currently the definition of more than 7000 triangle centers as well as the description of many properties 
of these points. Similarly, an encyclopedia of more than 700 cubics is maintained by Bernard Gibert and an 
encyclopedia of quadri-figures is maintained by Chris van Tienhoven.* 


l http://faculty.evansville.edu/ck6/encyclopedia/ETC.html. 
2 http://bernard. gibert.pagesperso-orange.fr/. 
3 http://www.chrisvantienhoven.nl/index.php/mathematics/encyclopedia. 


The website is available here: http://dpt-info.u-strasbg.fr/~narboux/CETC/about.html. 
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Table 1 Some ETC points 


Number Name Trilinear coordinates Barycentric coordinates 
Xı Incenter 1 a 
Xo Centroid l/a 1 
X3 Circumcenter cos(A) a? (a? — b? — c?) 
X4 Orthocenter sec(A) (a? +b? — c?) (a? — b? +c”) 
X5 Nine-point center cos(B — C) ab? — bt +a? + 2b*c? — cl 


The ETC consists in web pages containing hundreds of thousands of properties. But most of the results are pre- 
sented without proof nor reference. A large majority of the results have been generated by computer programs using 
symbolic and/or numeric manipulations. But the source code of the computations are not given in the encyclopedia 
so the reader can not reproduce the results and we can not verify if the result has been found by symbolic or numeric 
computations. How can one trust theses results ? What are the potential sources of errors ? 


1. There could be an error in the definition of the points in the Computer Algebra System (CAS).4 

2. There could be an error in the definition of the algorithms to check properties implemented in the language of 
the CAS. 

There could be a bug in the implementation of the CAS itself or in the theory behind the algorithms used. 

4. An error could have been introduced while copying the results from the computer software to the web-page. 


Our aim is to reduce these risks of errors. For this purpose we generate automatically the description of the points 
and we provide formal proofs. We created a certified version of ETC, such that each property of a point come along 
with a computer checked proof of its validity with the necessary non degeneracy assumptions. 

The work closest to ours is the work of Deko Dekov who uses a computer program to generate theorems in 
triangle geometry and present the results as a Computer Generated Encyclopedia of Euclidean Geometry. This 
work differs from ours because the proofs of the results are not given and the facts are not checked using a proof 
assistant. 

In the first section we provide an overview of the encyclopedia from the point of view of formal theorem proving. 
In the second section we describe our formalization using the Coq proof assistant of the required concepts. In the 
third section we explain how we obtained automatically thousands of properties and how we proved them. 


2 The Encyclopedia of Triangle Centers 


The triangle centers defined in the encyclopedia of Clark Kimberling are various. The first centers are well known 
points studied for centuries. The other points are defined by various means: most are defined by a traditional 
geometric constructions, some are defined as minimizer of some geometric quantity, there is even a point defined 
as a solution of an equation related to a physical property (X5626 is the center of maximal electrostatic potential 
inside a triangle ABC having a homogeneous surface charge distribution). Table 1 provides the list of the first 5 
points in ETC. a, b and c denote the side length of triangle ABC (a = |BC|, b = |AC| and c = |AB|). Figure 1 
depicts the first 4000 points of ETC. 

The points are defined by different means but are always manipulated using their homogeneous coordinates 
relative to the triangle ABC. The coordinates used are either the barycentric coordinates or the trilinear coordinates 
which are the relative distance from the point to the sides of the triangle ABC. These coordinates are expressed 


4 For example the definitions of points X559g to X 5602 contained an error in the version of ETC as of May 2015, the point was defined 
as a+ — a° (b + c)? — 4aS/r F 4R instead of a* — a? (b +c)? — AaSVrR + 4R2. 


5 http://www.ddekov.eu/j/contents.htm. 
6 http://www.ddekov.eu/el/ 
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Fig. 1 The first 4000 centers displayed (sketch generated using GeoGebra 5.0) 


A 


Fig. 2 Trilinear coordinates 


using the length of the sides: a, b, c or equivalently the measure of the three angles of triangle ABC. A given point 
P has trilinear coordinates any triple of the form ka’ : kb! : kc’ with k 4 0 and a’, b’ and c’ denote the distances 
from P to the side lines of ABC (Fig. 2). 

Kimberling defined a triangle center as a point X such that there is a function f of the side lengths of triangle 
ABC such that: X = f(a,b,c): f(b,c,a): f(c,a, b) and f should be homogeneous, there must be an integer n 
such that: Vabct,t > 0 => f(ta, tb, tc) = t” f(a, b, c) and symmetric in the sense that f(a, b,c) = f(a, c, b). 
The fact that f is a function of the side length ensure that the triangle centers are preserved by isometries (as the 
side lengths are preserved by isometries). The fact that f is homogeneous moreover ensures that triangle centers 
of similar triangles are preserved. 

To give to the reader an idea of what can be found in Clark Kimberling’s encyclopedia, we reproduce here the 
entry of the orthocenter (Fig. 3). For most points in the encyclopedia the trilinear and barycentric coordinates are 
given (often with several equivalent formulas). Then, some comments are given about the point in natural language. 
Finally, we have a list of properties. The lines that the point belongs to are listed. Then a list of relations between 
the points is given. For instance, X4 is the midpoint of the segment X3 X3g2. The amount of information available 
differs from one point to the other. 
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X(4) = ORTHOCENTER 


Trilinears sec A: sec B: sec C 
= cos A - sin B sin C : cos B- sin C sin A: cos C - sin A sinB 
= cos A - cos(B - C) : cos B - cos(C - A) : cos C - cos(A - B) 
= sin B sin C - cos(B - C) : sin C sin A - cos(C - A) : sin A sin B - cos(A - B) 


Barycentrics tan A : tan B : tan C 

Asa point on the Euler line, X(4) has Shinagawa coefficients (0, 1). 

X(4) is the point of concurrence of the altitudes of ABC. 

If you have The Geometer’s Sketchpad, you can view Orthocenter. 

X(4) and the vertices A,B,C comprise an orthocentric system, defined as a set of four points, one of which is the orthocenter of the triangle of the other three (so that each is the orthocenter of the other three 
Suppose P is not on a sideline of ABC. Let A'B'C' be the cevian triangle of P. Let D,E.F be the circles having diameters AA', BB',CC’. The radical center of D,E,F is X(4). (Floor van Lamoen, Hyacinthos #. 
Ross Honsberger, Episodes in Nineteenth and Twentieth Century Euclidean Geometry, Mathematical Association of America, 1995. Chapter 2: The Orthocenter. 

X(4) lies on the Thomson cubic, the Darboux cubic, the Napoleon cubic, the Neuberg cubic, and these lines: 

1,33 2.3 653 7.273 8.72 910 11,56 12,55 13,61 14,62 15,17 16,18 32,98 35,498 36,499 37,1841 39.232 42,1860 46,90 48.1881 49,156 51,185 52,68 54,184 57,84 64,1853 ı 
120,1292 121,1293 122,1294 123,1295 126,1296 127,1289 128,930 129,1303 130,1298 131,135 137,933 141,1350 145,149 147,148 150,152 155,254 162,270 165,1698 171,601 19. 


487,489 488.490 496.999 512,879 523,1552 542,576 569,1179 572.1474 574.1506 575.598 579.1713 580.1714 590,1151 608.1518 615.1152 616.627 617,628 653,1156 774,1254 8I 
1161,1163 1251,1832 1329,1376 1340,1348 1341,1349 1385,1538 1430,1468 1499.1550 17151730 1716,1721 1717,1718 1726,1782 3065.3464 3347,3472 3348.3355 


X(4) is the {X(3),X(5)}-harmonic conjugate of X(2). For a list of other harmonic conjugates of X(4), click Tables at the top of this page. 
X(4) = midpoint of X(i) and X(j) for these (ij): (3.382), (146.3448), (147,148). (149.153), (150.152) 


X(4) = reflection of X(i) in X(j) for these (i,j): (1,946), (2,381), (3,5), (5,546), (8,355), (20,3), (24,235), (25,1596), (40,10), (69,1352), (74,125), (98,115), (99,114), (100,119), (101,118), (102,124), (103,1 1¢ 
(1114,1313), (1.141.137), (1292,120), (1293121), (1294,122), (1295,123), (1296,126), (1297,127), (1298,130), (1299,135), (1300,136), (1303, 129), (1350,141), (1593,1595) 
X(4) = isogonal conjugate of X(3) 

X(4) = isotomic conjugate of X(69) 

X(4) = cyclocevian conjugate of X(2) 

X(4) = inverse-in-circumcircle of X(186) 

X(4) = inverse-in-nine-point-circle of X(403) 

X(4) = complement of X(20) 

X(4) = anticomplement of X(3) 

X(4) = complementary conjugate of X(2883) 

X(4) = anticomplementary conjugate of X(20) 

X(4) = eigencenter of cevian triangle of X(i) for I= 1, 88, 162 

X(4) = eigencenter of anticevian triangle of X(i) for 1 = 1, 44, 513 


X(4) = X(i) -Ceva conjugate of X(j) for these (i,j): 
(7,196), (27.19), (29,1), (92,281), (107,523), (264.2), (273.278), (275,6), (286.92), (393, 459) 


X(4) = cevapoint of X(i) and X(j) for these (i,j): 
(1,46), (2.193), (3.155), (5,52), (6.25), (11,513), (19,33), (30,113), (34.208), (37.209), (39,21 1), (51.53), (65.225), (114.511), (115,512), (116,514), (117,515), (118,516), (119,517), (120,518), (121,519), (1 


X(4) = X(i) -cross conjugate of X(j) for these (i,j): 
(3,254), (6.2), (19.278), (25,393), (33,281), (51,6), (52.24), (65,1), (113.403), (125,523), (185.3), (193.459), (225,158), (389.54), (397.17), (398.18), (407.225), (427.264), (512.112), (513,108), (523,107) 


X(4) = crosspoint of X(i) and X(j) for these (ij): (2.253). (7.189). (27.286), (92,273) 


X(4) = crosssum of X(i) and X(j) for these (i,j): 
(4,1075), (6,154), (25,1033), (48.212), (55,198), (56,1035), (71.228). (184,577), (185.417). (216.418) 


X(4) = crossdifference of any two points on line X(520)X(647) 


X(4) = X(i) -Hirst inverse of X(j) for these (ij): 
(1,243), (2,297), (3,350), (19,242), (21,425), (24.421), (25,419), (27.423), (28.422), (29.415), (420,427), (424.451), (459,460), (470,471), (1249.1503) 


X(4) = X(i) -aleph conjugate of X(j) for these (ij): (1.1047). (29.4) 


X(4) = X(i) -beth conjugate of X(j) for these (i,j): 
(4.34), (8.40), (29.4), (162.56), (318.318), (811.331) 


Fig. 3. An example of one of the 7000 entries in Clark Kimberling’s encyclopedia: point X4 orthocenter 


2.1 Categories of Points 


We can distinguish triangle centers depending on how they can be expressed using trilinear coordinates, this will 
be important to obtain proofs. 

About 87 % of the points can be defined by trilinear coordinates which are rational fractions in a, b and c. The 
definition of quite a few points (about 12 %) involve square roots, either square roots of constants or also square 
roots of rational fractions in a, b and c. There are only a few occurrences of nested square roots in the original list of 
coordinates collected by Peter Moses. One definition involve a constant (v 25 — 104/5 in the definition of X 1139) 
but the other occurrences can be simplified to suppress the nested radicals. Some points (for instance X506 and 
X507) contain cube roots or fourth roots (.!/3, .!/4) The definition of about 1 % of the points contains trigonometric 
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Table 2 Conway’s notations 


Symbol Description Definition 
À Tangle arei V(iatb—c\a—b+c)\(-a+b+cy\a+b+c) 
4 
a? +b? +c? 
SW LOO 
2 
b+c-—-a 
sa — 
2 
b c+a-—b 
s ae 
2 
a+b—c 
sc — 
2 
; tote 
s Half-perimeter Gr 
J Ja® — ath? — a? bt + b6 — atc? + 3a2b2c? — btc? — a?c — b?ct + cô 
abc 
af — a?b? + bt — a?c? — bc? + ct 
e 
a2b2 + a?c? + b2c2 
b ~ 
R Radius of the circumcircle = 
V(iatb+cy(b+c—a\(c+a—by(a+b—c) 
1 (b+c : b tb 
r Radius of the incircle L caleta Xa c) 
2 a+b+c 


functions which can not be eliminated : sin, cos, tan, sec, csc, sec, cot. For many points (such as the circumcenter 
or orthocenter) the definition using trigonometric functions is short but they can be eliminated to ease computations 
using the law of cosines (cos (A) = (b? + c? — a*)/(2bc)). 

The definition of some points make use of the value of the angles (A, B, C) without using any trigonometric 
function. For example X360, Hofstadter zero point, has trilinear coordinates: A/a : B/b : C/c. For very few points 
(X368, X369, X370, X 1144, X3232, X5373, X5626) we do not have the explicit coordinates, thus we do not treat those 
points. In order to obtain smaller expressions for the coordinates, some notations (known as Conway’s notations) 
are used, Table 2 lists those notations. 


3 Formalization of the Main Concepts Using Coq 


In this section, we describe the formalization in Coq of the concepts we need. The results given in this section are 
stated and proved interactively using the Coq proof assistant (even if we use some automatic tactics). In the next 
sections, we will describe the facts which have been found and proved completely automatically. 

For the Coq formalization we start with the type of the points, we use a record with three real numbers to represent 
the homogeneous coordinates. We don’t specify in the type the exact kind of coordinate we use, it can be either 
trilinear or barycentric coordinates. This type will also be used to represent lines because of the duality between 
points and lines using homogeneous coordinates. 


Record Point : Set := 


cTriple 
{ 
X Re 
x R; 
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Note that a unique point can be represented using different homogeneous coordinates, hence we define point 
equality as: 


Definition eq_points P Q := X(Q)*Y(P) = X(P)*Y(Q) /\ X(Q)*Z(P) = X(P)*Z(Q). 


We also define a type class representing the given base triangle defined by three real numbers: the lengths of the 
sides a, b and c. We assume that the sides are strictly positive and verify the triangle inequality. 


Class triangle_theory := { 
a:R; 

b:R; 

GR; 
apos: a > 
bpos: b > 
cpos: c > 
ineq_1 
ineq_2 
ineq_3 : -a + 


In this setting, the three reference points can be defined as:’ 


Definition pointA := cTriple 1 0 0 
Definition pointB := cTriple 0 1 0. 
Definition pointC := cTriple 0 0 1 


As many theorems are incorrect for points on the side lines of triangle ABC, to simplify our statements we 
introduce a predicate: 


Definition is_not_on_sidelines P := 
match P with 

cTriple X Y Z => X<>0 /\ Y<>0 /\ Z<>0 
end. 


To ease the definition of triangle centers we define a constructor which take as an argument an homogeneous 
function h and produce the associated triangle center: 


Context ‘{M:triangle_theory}. (* We assume we have a triangle theory. *) 


Definition cPointh h := cTriple (ha b c) (hbca) (hca b). 


Our development uses trilinear coordinates, but the definition of the coordinates of the points that Peter Moses 
kindly gave us, uses barycentric coordinates. Fortunately it is easy to convert from one to the other. For more 
convenience we define an alternative constructor to obtain trilinear coordinates from a function for barycentric 
coordinates: 


Definition cPointhb h := cTriple ((h a b c)*b*c) ((h b c a)*a*c) ((h c a b)*ax*b). 


Using these constructors we define 6000 points, we list here the first five points: 


Definition X_1 := let h_x_1l abc := a in cPointhb h_x_l. 

Definition K_2 := let h_x_2 abe := 1 in cPointhb h x 2. 

Definition X_3 := let h_x_3 a b c := a^2*(a^2-b^2-c^2) in cPointhb h_x_3. 

Definition X_4 := let hx 4a bc := (a*2+b*2-c%2)*(a*2-b*2+c*2) in cPointhb h_x_4. 
Definition X_5 := let h_x_5 a b c := a^2*b^2-b^°4+a^2*c^2+2*b^2*c^2-c^4 in cPointhb h_x 5 


7 We name the points like that to keep the letters A B C for other purposes. 
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3.1 Incidence Geometry 


A very important predicate in geometry is the collinearity predicate. Using homogeneous coordinates three points 
are collinear if their determinant is zero, we can define the col predicate as: 


Definition col ABC := 
X(A) *Y¥(B)*Z(C) + Y¥(A)*Z(B)*X(C) + Z(A) *X(B) *Y(C) 
= Z(A)*Y(B)*X(C) = Y(A)*X(B)*Z(C) = X(A)*Z(B)*Y(C) = 0. 


The line going through P and Q can be defined as a triple (the three coefficients of the equation of the line): 


Definition line P Q := 
match P,Q with 
cTriple xl yl z1, cTriple x2 y2 z2 => 
cTriple (y1*z2-y2*z1) (z1*x2-z2*x1) (x1l*y2-x2*y1) 
end. 


Using the duality, the point of intersection of two lines can use the same definition, but we provide an alternative 
definition to increase readability: 


Definition inter L M := line L M. 


3.2 Euclidean Geometry 


X1 y1 Z1 
Two lines of coordinates (x1, y1, Z1) and (x2, y2, z2) are parallel if |x2 y2 z2| = 0: 
abc 


Definition parallel L M := 
match L,M with 
cTriple xi yl z1, cTriple x2 y2 z2 => 
x1l*y2*c+x2*b*zlta*yl*z2 -a*y2*z1-x1*b*z2-x2*yl*c = 0 
end. 


We can then easily prove the triangle midpoint theorem to check our definitions: 


Lemma triangle_midpoint : forall ABC, 
parallel (line (midpoint A B) (midpoint A C)) (line B C). 


Now, we can define a predicate for perpendicular lines. We adapt the formula which can be found in the litera- 
ture [6] to get rid of trigonometric functions and divisions: 


Definition perp LM := 
match L,M with 
cTriple x1 yl z1, cTriple x2 y2 z2 => 
2*a*b*c* (x1*x2 + yl*y2 + z1*z2) 
—(yl*z2 + y2*z1)* a* (b*%2+c%2-a%2) 
—(z1*x2 + zZ2*x1)* b* (c*2+a%2-b*2) 
—(xl*y2 + x2*y1)* c* (a%2+b*2-c%2) = 0 
end. 
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Fig. 4 Cevian triangle and pedal triangle 


Then, we can check that the circumcenter (X3 in Kimberling’s numbering) belongs to the perpendicular bisectors: 


Lemma circumcenter_is_on_perp_bisect 
perp (line X_3 (midpoint pointA pointB)) (line pointA pointB). 


We define the foot of the altitudes: 


Definition T_A := inter (line pointB pointC) (line pointA X_4). 
Definition T_B := inter (line pointA pointC) (line pointB X_4). 
Definition T_C := inter (line pointA pointB) (line pointC X_4). 


With these definitions we can check that the orthocenter is properly defined: 


Lemma orthocenter_check : 
perp (line X_4 T_A) (line pointB pointcC) /\ 
perp (line X_4 T_B) (line pointA pointcC) /\ 
perp (line X_4 T_C) (line pointA pointB). 


3.3 About Triangles 
3.3.1 Cevian and Anti-Cevian Triangles 


The concept of Cevian triangles will be useful in the rest of this paper. In this paragraph we define the Cevian and 
anti-Cevian triangles of a point P. We adopt a definition based on coordinates and we will prove afterward that 
the definition is correct. Recall the Cevian triangle of a point P in a triangle ABC is the triangle formed by the 
intersections of the three lines joining the vertexes of ABC to P and the sides of ABC (Fig. 4). A’B’C’ is the 
anti-Cevian triangle of P in ABC, if ABC is the Cevian-triangle of P in ABC. 
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Definition cevian_triangle P := 

match P with cTriple p q r => 

(cTriple 0 q r, cTriple p 0 r, cTriple p q 0) 
end. 


Definition anti_cevian_triangle P := 
match P with cTriple p q r => 
(cTriple (-p) q r, cTriple p (-q) r, cTriple p q (-r)) 
end. 


But these definitions allow us to define only the Cevian triangle of P w.r.t. the base triangle. If we want to 
consider the Cevian-triangle of P in any triangle ABC, we need a more general definition: 


Definition cevian_triangle_gen P ABC := 
( (inter (line A P) (line B C)), 
(inter (line B P) (line A C)), 
(inter (line C P) (line A B))). 


We can convince ourselves that particular case for the base triangle ABC is consistent with the general definition: 


Lemma cevian_triangle_ok : forall P, 
eq triangles (cevian_triangle P) (cevian_triangle_gen P pointA pointB pointcC). 


We can prove also that we really defined the Cevian triangle: 


Lemma cevian_triangle_gen_ok 
forall PA BC, 
match cevian_triangle_gen P A B C with 
(X,Y,Z) => col A P X /\ col BP Y /\ col CPZ /\ 
col X BC /\ col YAC /\ col ZAB 
end. 


Finally, we can name some well known triangles: 


Definition incentral_triangle := cevian_triangle X_1. 
Definition medial_triangle := cevian_triangle X_2. 
Definition orthic_triangle := cevian_triangle X_4. 
Definition intouch_triangle := cevian_triangle X_7. 
Definition extouch_triangle := cevian_triangle X_8. 


3.3.2 Pedal Triangle 


The Pedal triangle A p BpCp of a point P is the triangle formed by the projections of P on the side lines (Fig. 4), 
we can define its coordinates: 


Definition pedal_triangle P := 
match P with cTriple p q r => 


(cTriple 0 (qt+tp* (-a*2-b*2+c%2)/(-2*a*b)) (r + p*(-a^2+b^2-c^2)/(-2*a*c)), 
cTriple (p + q*(-a^2-b^°2+c^2)/(-2*a*b)) 0 (r+q*(a^2-b^2-c^2)/(-2*b*c)), 
cTriple (p + r*(-a%2+b*2-c%2)/(-2*a*c)) (qtr* (a*2-b*2-c%2)/(-2*b*c)) 0) 


end. 
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Fig. 5 Perspective triangles 


As previously, we prove afterward that our definition make sense: 


Lemma pedal_triangle_ok : forall P, 
match pedal_triangle P with 
(A1,B1,C1) => perp (line P A1) (line pointB pointCc) /\ 
perp (line P B1) (line pointA pointc) /\ 
perp (line P C1) (line pointA pointB) 
end. 


Obviously, the Cevian triangle of the orthocenter is the same as the pedal triangle: 


Lemma cevian_triangle_of_orthocenter_is_pedal 
eq_triangles (cevian_triangle X_4) (pedal_triangle X_4). 


Moreover, the medial triangle is the pedal triangle of the circumcenter: 


Lemma pedal_triangle_of_circumcenter_is_medial 
eq_triangles medial_triangle (pedal_triangle X_3). 


The intouch triangle is the pedal triangle of the incenter: 


Lemma pedal_triangle_of_incenter_is_intouch 
eq_triangles intouch_triangle (pedal_triangle X_1). 


3.3.3 Perspective Triangles 


Two triangles ABC and A’B’C’ are said to be perspective from O if their three pairs of corresponding vertexes are 
joined by lines which meet in O (Fig. 5): 


Definition is_perspector O T1 T2 := 
match T1,T2 with 
(A,B,C), (A , B,C") = 
col OA A’ /\ col OBB’ /\ col OCC’ 
end. 
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Table 3 The transformations formalized 
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Name#Arity Definition 
—a b c m 

Anticomplement#1 TUPE Oger ap—bg+er aptbg—er 

a 
Ceva_conjugate#2 u(—qru +rpv + pqw) v(qru —rpv + pqw) w(qru +rpv — pqw) 

bq + cr ap + cr ap + bq 

Complement#1 — — 

a b c 

u v w 


Cross_conjugate#2 


Cross_difference#2 
Cross_sum#2 
Cross_point#2 


Daleth_conjugate#2 


Hirst_inverse#2 
Iso_conjugate#2 


Isogonal_conjugate#! 


Isotomic_conjugate#| 
Line_conjugate#2 


Midpoint#2 


Reflection#2 


Waw_conjugate#2 


—pvw+qwu+ruv 


qu—ruv 
qu +rv 
pu(ru +qw) 


let h(p,q,r,u, v, w) = pÈ 
h(p, q,r, u,v, w) 


qru? — vwp? 
qrvw 


letk = au + bv 
kp + hu 

(ap — bq — cr)u 
+2p(bv + cw) 


pvw —qwu+ruv 
ru — pw 
ru + pw 
qu(pw +ru) 
2 2u v w. 
y+ ul ) in 
q r 
h(q,r, p, v, w, u) 
rpv? — wuq 
rpwu 
1 


q 
1 
qb? 
q(w? +u?) 
—v(rw + pu) 


w 


cr in 
cw in 
kq + hv 
(bq — ap —cr)v 


+2q (au + cw) 
2724 


2p°(ry — qz)? — pqr’ xy 


pvw +qwu — ruv 
pv— qu 
pv + qu 
rw(qu + pv) 


h(r, p,q, w, u, v) 
pqw? — uvr? 
pquv 


r(u2 + v?) 
—w(pu + qv) 


kr + hw 
(cr — ap — bq)w 
+2r (au + bv) 


let h(p, q,r, x,y,z) = p(x?q 


h(q,r, p, v, w, u) 
2 2 


h(p,q,1,u, v, w) 


Zayin_conjugate #2 


h(p, q,r, u, v, w) h(q,r, p, v, w, u) 


For example, P is the perspector of any triangle ABC and its Cevian triangle in ABC: 


let A(p, q,r, x, y,7) = PO +2)? — ry? — qz? + (p —r)xy + (p — q)xz in 


h(r, p,q, w, u, v) 


h(r, p,q, w, u, v) 


Lemma cevian_triangle_perspector : 
forall A B C P, is_perspector P (cevian_triangle_gen P A B C) 


3.4 The Geometric Transformations 


We define the geometric transformations using homogeneous coordinates. Table 3 lists all the transformations we 
can deal with. The definitions of these transformations are taken from Kimberling’s glossary or Paul Yiu’s book 
[15]. We could have given the definitions using geometric constructions, but we adopt algebraic definitions for two 
reasons. First, with algebraic definition, the proofs are relatively easy to obtain. Second, the geometric definitions 
often require a proof to certify the definition make sense. For example, the geometric definition of the centroid 
requires to prove that the medians intersect. Still, we want to make sure we did not make any mistake while giving 


our definitions. To check our definitions, we prove the geometric properties. 


Ceva Conjugate. As a first example, let’s consider the P-Ceva conjugate of Q. By definition it is the perspector of 


the Cevian triangle of P and the anti-Cevian triangle of Q but we define it as: 


68 J. Narboux, D. Braun 


Definition ceva_conjugate P U := 

match P,U with 

cTriple p q r, cTriple u v w => 

cTriple (u*(-q*r*u+tr*p*v+p*žq*w)) (v* (q*r*u-r*p*vtp*q*w)) (w* (q*r*utr*p*v-p*q*w) ) 
end. 


To prove that our definition is consistent with the geometric definition, we can prove that the P-Ceva conjugate 
of Q is the perspector of the Cevian triangle of P and anti-Cevian triangle of Q (Fig. 6): 


Lemma ceva_conjugate_property : 
forall P Q, is_perspector (ceva_conjugate P Q) (cevian_triangle P) (anti_cevian_triangle 
Q). 


Isotomic Conjugate. A second example is the isotomic conjugate. Given a point P, let A; B,C, be the Cevian 
triangle of P in ABC. Let A2 (resp. B2, C2) be reflection of A; (resp. B1,C1) in the midpoint of segment BC (resp. 
AC,AB). The lines AAz, B Bz and CC? meet at a point P’ called the isotomic conjugate of P (Fig. 7). We define 
it using coordinates as: 


Fig. 6 The Ceva conjugate 


Fig. 7 The isotomic conjugate 
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Definition isotomic_conjugate P := 

match P with 
cTriple p q r => 
cTriple (1/(p*a*2)) (1/(q*b*2)) (1/(r*c%2)) 
end. 


Then to prove that our definition is correct, we express the fact that if we reflect the Cevian triangle of P in the 
midpoints of the sides of ABC we obtain the Cevian triangle of its isotomic conjugate. Note that we need to assume 
that P is not on any sideline of triangle ABC. 


Lemma isotomic_conjugate_property 
forall P, 
is_not_on_sidelines P -> 
match cevian_triangle P with 
(A1,B1,C1) => 
match cevian_triangle (isotomic_conjugate P) with 
(A2,B2,C2) => eq points (midpoint pointA pointB) (midpoint Cl C2) 
end 


Easily we can show also that this transformation is an involution: 


Lemma isotomic_involution : forall P, 
is_not_on_sidelines P -> 
eq_ points (isotomic_conjugate (isotomic_conjugate P)) P. 


Complement and Anti-Complement. U is the complement of a point P if it belongs to line X2P (where X2 is the 
centroid) and X2 trisects segment PU and X2 is closer to P than U. We define it as: 


Definition complement P := 

match P with 
cTriple p q r => 
cTriple ((b*q+tc*r)/a) ((a*pt+tc*r)/b) ((a*ptb*q) /c) 
end. 


Definition anticomplement P := 

match P with 
cTriple p qr => 
cTriple ((-a*p+b*qtc*r)/a) ((a*p-b*qt+tc*r)/b) ((a*p+b*q-c*r) /c) 
end. 


We can prove that they are inverse function of each-other: 


Lemma anticomplement_complement_property 
forall P, eq_points (anticomplement (complement P)) P. 


and that the centroid belongs to the line: 


Lemma complement_property 
forall P, col P (complement P) X_2. 


Reflection and Midpoint. Following the formulas given in Table 3, we can define the midpoint of a segment and the 
reflection of point wrt. another point, then we can check that: 


Lemma reflection_midpoint : forall A B, 
eq_ points (midpoint A (reflection A B)) B. 
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As an example, we can also prove that the center of the nine point circle is the midpoint of the segment X3X4 
(X3 is the circumcenter and X4 the orthocenter): 


Lemma nine_point_circle_bisects_OH: 
eq points X_5 (midpoint X_3 X_4). 


4 Finding Conjectures and Automating their Proofs 


In the previous section, we listed our definitions and presented some results that we proved interactively just to 
verify our definitions. In this section we describe the algorithms we used to find and prove automatically some 
properties about triangle centers. 


4.1 Finding Properties 


There are so many properties presented in Clark Kimberling’s encyclopedia that obtaining an exhaustive list of 
these properties is difficult. Indeed, it would require to parse the information presented on the web pages, but the 
computer generated pieces of information are inter-mixed with human generated comments in such a way that it 
would make this parsing phase non-trivial. Hence, we chose not to check the properties displayed on the page, but 
to generate ourself these properties. Fortunately most of the properties displayed in the encyclopedia are relatively 
easy to obtain by computer means. In this section, we describe how we obtained these properties. 

We address two main kind of problems: 


1. Find all the collinearity properties in the list of points. Namely find all the triples (7, j,k) such that 
Col(X;, Xj, Xx). 

2. Given a function f which takes n points as arguments and return a point, find all n + 1-tuples (i, ji, ..-, jn) 
such that X; = f(X js- Xn) 

We want to find properties which are valid in any triangle, but to simplify the computations we consider one 
example triangle. If the property does not hold in this triangle, it does not hold in general, if it holds in this 
particular triangle then there are good chance that it holds in general. In practice the coincidences are very rare 
[but Francisco Javier recently pointed that points X3g75 and X7292 have exactly the same coordinate in the triangle 
(a, b, c) = (6, 9, 13)]. 


4.1.1 Finding Collinearity Properties 


Letn be the number of points in ETC. We want to find all the triples (i, j, k) such that Col(X;, X j, Xk). A naïve algo- 
rithm would require O(n*) operations. We implemented a better algorithm by using an associative array T to store 
for each line the set of points belonging to this line. For all every couple of points (X;, X j) we compute the equation 
e of the line X; X; and we add the points X; and X; to T[e]. This algorithm has a complexity O(n?) on average. 

The main difficulty concerns the numerical computations. Our first prototype implemented using Maple™was 
producing both false positive and false negative: some points appeared to be collinear whereas there were not 
(even in the special case of the chosen triangle) and also some points which are collinear were reported not to 
be collinear. This is due to the fact that Maple™(or other similar software) does not compute itself the precision 
needed to obtain a significant result. In a second prototype, we use what we call a symbolic-numeric approach. 
Instead of representing the coordinates of the points using big-floats approximations, we used symbolic expressions 
involving constants, arithmetic operations (+,—,*,/), square roots and trigonometric functions. Deciding equality of 
two expressions involving radicals is possible but not trivial: for example V3 + 2/2 = 1 — V2. To decide such 
an equality involving only constants we use the radnormal function of Maple™. But we are not aware of any 
algorithm dealing with the trigonometric functions. To prove the facts, we need an algorithm to normalize sym- 
bolic expressions containing radicals, an algorithm has been proposed by Rybowicz which generalize an algorithm 
introduced by Zippel [13,16] but we are not aware of any implementation in a CAS. 
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4.1.2 Finding Properties About Transformations 


To find the properties about transformations we use a brute-force algorithm. Recall that, given a function f 
which takes n points as arguments and return a point, we need to find all n + 1-tuples (i, ji, ..., Jn) such that 
Xi = f(Xj,,..., X jp). For every tuple of points (X;,,..., Xj) we compute X = f(X;,,..., X ja) and we check 
if X belongs to the (sorted) list of points. The complexity of this algorithm is O(k"/n(k)) where k is the number of 
points in the encyclopedia. 


4.2 Proving the Properties 


To prove the different properties in Cog, we make heavy use of the ring and field tactics of Coq first introduced by 
Delahaye and Mayero and improved by Grégoire and Mahboubi [5,7]. The idea of these reflexive tactics is the follow- 
ing: instead of convincing Coq that two expressions representing the same polynomial are equal by repeated applica- 
tions of commutativity, associativity and distributivity properties, an algorithm to compute a normal form is written in 
Coq’s language and proved correct. For more information about how to develop reflexive tactics see Chapter 16 of [1]. 

Proving properties involving only points whose coordinates are rational fractions is straightforward, we just 
need to unfold all the definitions to obtain a goal which can be solved by the field tactic. We do not treat the points 
whose coordinates involve trigonometric functions. For the points with square roots in their coordinates we use the 
following heuristics: 


1. We try to simplify square roots by applying the following lemma repeatedly (where S n stands for n + 1): 


Lemma sqrt_pow : forall x n, 0 <= x -> 
(sqrt x)*(S (S n)) = x* (sqrt x)“n. 


2. We normalize all notations before unfolding them. For each notation, we have a lemma stating a permutation 
property, such as: 


Lemma Delta_perm_1 : DeltaMaj a b c = DeltaMaj b c a. 


3. We normalize all square roots of rational fractions appearing in the expression. We have a tactic 
normalize_all_sqrt which computes the normal form of every square root appearing in the goal. Here 
is an example of its use: 


Goal sqrt ((at+b)*2) + sqrt (a+b) = sqrt (bta)+sqrt(a*2+b*2+2*a*b). 
Proof. 

normalize_all_sqrt. 

ring. 

Qed. 


Non Degeneracy Conditions. The non degeneracy conditions needed for proving the properties are that all the 
denominators appearing in the statements are different from zero and the the argument of any square root is posi- 
tive, they can be generated automatically. 


4.3 Results 


Figure 8 shows a snapshot our website displaying the pieces of information about the point X 1. Each fact is displayed 
in three colors: 
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X(1) 


gare Bs d 


X(1) is isogonal conjugate of X(1) (EIS 

X(1) is isotomic conjugate of X(75) 

X(1) is cyclocevian conjugate of X(1029) (aay 

X(1) is complement of X(8) EE 

X(1) is Ceva-conjugate of (i,j 

e X(1) is cross-conjugate of 

X(1) is cross-point of (i,j) for Și 6 7 (7,369) (7,370) (7,1144) (7,1171) (7 7, 1525) { (7, 1564) (7, 3130) (7 

(7, 385 7, 5373) (7, 5394) 0 J (8,280) (10,502) (21,29) (59,110) (75,92) (81,86) (82,3112) (88,679) (99,1016) (100, 765) 


ec eee 


(1222, 1476) (4567, 4596) (4618, 5376) 


57) (4, 1148) (6,55) (9, 3158) 
(50, 215) (56,221) (57,1419) (58,5 7 | 
(492, 1564) (192.3130) { ( 50 5373 5: 192, 5403 44, 513) (354 (385, 43 
(512, 1015) (614, 4319) (647. | 2310) | : | j 2 (774,820) (798, 4117) 
(810, 3708) (984, 3795) (999 | 1 8) { 9 2 ( 238 69, 3725) (2594, 2599) 


. 


(1,650) (1.652) (1.654) (1,656) (1,657) (1,659) (1,661) (1,672) 
(1,910) (1,1155) (1, 1491) 5 (1,1755) (1,21 2 (1,2183) ( 2 ( ye 
(1, 2236) (1 j 2244) 1 ) 2252) 


2484) (1,2509) (1,2511) (1,2515) 


( 2631) (1, 2637) (1.2641) (1, 2642) 
(1, 3013) (1, 3287) (1,3330) (1, 3768) (1.4394) (1, 4724) 1 2) 4784) (1 ‘90 4893) 4979) (2,3008) (7,516) (75,1733) (82, 3405) 
(679, 1) (765.1) (1910, 3404) 


X(1) is Waw conjugate of (i,j) for 
X(1) is line conjugate of (i,j) for 580 pairs: 


Fig. 8 A snapshot of the first entry of the Certified ETC 


Table 4 Statistics 


Name Number of facts proved Number of conjectures Percentage 


ceva_conjugate 7592 13091 58 
Complement 906 932 97 
Cross_conjugate 2249 5483 41 
Cross_sum 4989 8552 58 
Cross_point 3381 8052 41 
Daleth_conjugate 8939 16216 55 
Hirst_inverse 9130 16383 56 
Isogonal_conjugate 1905 1954 97 
Isotomic_conjugate 604 616 98 
Line_conjugate 40136 57074 70 
Waw_conjugate 83 805 10 
Total 79,914 129,158 62 


Green We have a proof verified by the Coq proof assistant of this fact (assuming the non degeneracy conditions). 
Orange We have a Maple™ script to check this fact. 
Red This result is only a conjecture obtained by numerical computations. 


Table 4 shows the number of conjectures and the number of properties proved. It takes about 20 computer.days for 
Coq to compile these files. 


5 Conclusion 


We presented our first attempts toward a certified version of the encyclopedia of triangle centers. This experiment 
shows that the approach based on trilinear coordinates is efficient and can be used within a proof assistant. Compared 
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to other automated deduction approaches such as Wu’s or Grobner’s bases methods [4,8,14], this method consists 
in keeping always the explicit coordinates of the points, in some sense the triangulation of the algebraic system 
is computed in advance. This saves a lot of computations in this context where we want to prove many properties 
about the same given points. This work also demonstrates that the implementation of the Coq proof assistant scales, 
and can manage thousands of definitions. 

This work can be extended in several directions. 

First, it would be interesting to connect our proofs about trilinear coordinates to our development in the context 
of the geometry of Tarski [2,3,12]. For example we could prove that the definition of the collinearity predicate 
given in this paper is equivalent with the geometric definition given in the context of Tarski’s axiom system. 

Second, we could study how to find and prove automatically more properties. For example we could enhance our 
algorithm to find properties of some well known triangles, such as ‘the incenter is the circumcenter of the intouch 
triangle’. 

Third, we could improve the user interface. It would be interesting to have an interactive website allowing the 
user to ask questions in a simple syntax, to find automatically if a theorem is already known and to propose additions 
to the encyclopedia. 

Fourth, we could generalize this work and try to certify the results presented in the encyclopedia of conics and 
encyclopedia of quadri-figures. 
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